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Abstract 

We study single field inflation in noncommutative spacetime and compute two-point and three- 
point correlation functions for the curvature perturbation. We find that both power spectrum 
and bispectrum for comoving curvature perturbation are statistically anisotropic and the bispec- 
trum is also modified by a phase factor depending upon the noncommutative parameters. The 
non-linearity parameter fj^i is small for small statistical anisotropic corrections to the bispec- 
trum coming from the noncommutative geometry. There is a scale dependence of fisiL due to the 
noncommutative spacetime which is different from the standard single field inflation models and 
statistically anisotropic vector field inflation models. The anisotropy of the two-point function and 
the running of the f^L parameter will be tested precisely from the data which will be released soon 
from the PLANCK experiment. 

PACS numbers: 98.80.Cq, ll.lO.Nx, 98.70.Vc 
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1. INTRODUCTION 



Inflation (H not only solves the various puzzles of the Big-Bang theory, but it also provides 
seeds for the temperature anisotropy of the cosmic microwave background (CMB) radiation 
and structures in the universe. In the standard inflationary scenario, the potential energy of 
a scalar fleld called " inflaton" dominates the energy density of the universe and its quantum 
fluctuations generate perturbations in the metric causing small inhomogeneities in the early 
universe which give rise to CMB anisotropy and structures in the universe. 

Inflation predicts nearly scale invariant, adiabatic and gaussian perturbations. The flrst 
two are in excellent agreement with the observations of CMB anisotropy and polarization 
by COBE j2[, WMAP and other ground based and satellite based experiments, but 
the test of gaussian statistics of the perturbations is controversial and is the major goal of 
ongoing and future observations like Planck [6], CMBPOL [3] and Euclid satellite [8|. 

The non-gaussianity in CMB can be primordial or can be generated due to secondary 
sources (see j9| for detailed review). The primordial non-gaussianity arises due to the inter- 
action terms in the scalar potential and non-linearities of the gravity, where the latter effect 
is dominant than the former. The magnitude of the non-gaussianity in standard single-fleld 



inflation comes out to be small and of the order of slow-roll parameters [10 . 

Inflation occurs at a very high energy and it stretches out very small scales, of the order of 
Planck length, to the current hubble scale due to superluminal expansion of the universe. So. 
it provides window to see the new physics at the Planck scale at which quantum corrections 
to the gravity becomes important. These new effects can give rise to large primordial non- 
gaussianity that can be seen by Planck. 

Spacetime noncommutativity (see for review) is one of such modiflcations at high- 
energy, which is well motivated by quantum gravity and string theory. Modiflcations to 
the power spectrum of scalar perturbations during inflation and its effects on CMB due 
to noncommutative geometry has been studied in many places [l2-l^|. In this paper we 
compute the three-point correlation functions of the curvature perturbation and hence the 
non-linearity parameter fj^^ determining the primordial non-gaussianity using the noncom- 
mutative quantum fleld theories related to deformed Poincare symmetry [l2|. In this ap- 
proach of noncommutative geometry quantum flelds follow twisted statistics, as implied by 
the deformed Poincare symmetry in quantum theories. Non-gaussianity due to noncommu- 
tative spacetime of [l2[ has also been calculated in 16] using the 6N formalism which is 
used to calculate the local non-gaussianity and treats the comoving curvature perturbations 
as classical. But, here we compute the two-point and three-point function using Maldacena 
approach [lo| which is based on the second order perturbation theory and is valid for all 
kinds of non-gaussian effects. 



2. QUANTUM FIELDS IN NONCOMMUTATIVE SPACETIME 

At the energy scale of inflation, the noncommutativity of spacetime, which is motivated 
by Heisenberg uncertainty principle and Einstein's general relativity, can play a crucial role. 



The spacetime noncommutativity can be represented by the commutation relations [H 



[x^, X/y^ i9^u (1) 
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where 9fj_,y is a real antisymmetric matrix with constant elements and are the coordinate 
functions of the chosen coordinate system: 

~ ^^i- (2) 

The relation ([1]) holds only in special coordinate systems and looks quite complicated in 
other coordinates. The natural choice of the coordinate for cosmological applications is the 
comoving frame, where the galaxies are freely falling. We chose 9f^,y to be constant in the 
comoving coordinates as it makes our computation convenient and the physical scales can 
be computed by using the relations 

= am., e^{t) = a'm, (3) 

Our choice is opposite to the choice of [13|] where they assume 6^'^ to be constant in time. 

Due to spacetime noncommutativity, the usual quantum fields are deformed and can be 
given in terms of undeformed quantum fields as [ll| 

cPe = 0oe^^^^ (4) 

where d A P = h^B^'^P^, and P,^ represents the field momentum operator. The product of 
the deformed (twisted) quantum fields at the same spacetime point is represented by the 
star product given as 



x=y 

In the following sections we will make use of these relations to calculate two and three-point 
correlation functions of the comoving curvature perturbations. 



3. TWO-POINT AND THREE-POINT FUNCTION WITH NONCOMMUTATIVE 
SPACETIME 

3.1. Background 

The action of a single scalar field minimally coupled with gravity is 

S = l d'x^g(^^R + Cy (6) 
Here R is the Ricci scalar and C is the Lagrangian for the scalar field i.e. 

C = ~g^''d,cpd.c^~V {<!>). (7) 

Noncommutativity doesn't change the classical background so all the background dynamics 
will be similar to the standard case. We take the metric signature (— , +, +, +) and work in 
the units where Mp = 1. The background geometry of the homogeneous isotropic universe 
is described by the FRW metric 

ds^ = -dt^ + a?{t) [dx^ + dy"^ + dz^) . (8) 
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For a scalar field dominated universe the Friedmann equations are given as 



p + 3H{p + p) =0 

H = -^{p + p) (9) 

The equation of motion for the scalar field is given as 

4> + 3H^+^ = 0. (10) 
acp 

During inflation, the potential energy of the scalar field dominates the total energy density of 
the universe and the dynamics of the scalar field is governed by slow-roll parameters defined 
as as 

H 

V = — (11) 

Here we follow the definition of r] as in [l8| , which is different from the definition using 
scalar field potential {rjy = -^) and rj = —2r]y + 4e. 

3.2. Perturbations and ADM formalism 

The quantum fluctuations in the scalar field 50(x, t) generated during inflation are coupled 
to the perturbations in the metric through Einstein's equation. Inflation gives rise to scalar 
and tensor perturbations in the metric and the scalar part is written as 

ds^ = - (1 + 2$) dt^ + 2a^{t)B^idx'dt + a^{t) ((1 - 2^) % + 2E^ij) dx^dx^ (12) 

Here we have four scalar degrees of freedom in the metric and one in the scalar field which 
can be reduced to three by using gauge transformations. We can again use the constraint 
equations derived from the perturbed Einstein's equation and describe the scalar perturba- 



tions in terms of the curvature perturbation defined as [19 



C = -VP - ^50. (13) 

This variable is gauge invariant and is conserved on super-horizon scales. One can write the 
action (E]) in terms of C, and it turns out quadratic in C,. To do the perturbation theory in 
higher order it is convenient to use ADM formalism where the metric can be written as jl^ 

ds^ = -N^df + hij {dx' + N'dt) [dx^ + N^dt) (14) 
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where is laps function, Ni, Nj are shift vectors and hij is metric of three-dimensional 
hypersurface of constant time. Here and Ni appear as Lagrangian multipliers in the 
action so one can solve there constraint equation of motion and substitute the solution back 
into the action. This simplifies the tedious calculations needed while working with ( |T2|) . 
Now we chose comoving gauge Sep = to do our calculation and in this gauge we can use 



non-linear generalization of ( [21[ and define the gauge as [10|, 118 



(15) 



hij = a^e^^5ij, 5(j) = 
With this gauge choice the action ([6]) with the metric ( IT4l) becomes 

S=^J dtd^x^/h (^iVi?(3) - 2iVy (0) + N-^^"^ + N'^ {EijE'^ - E' 

Here R^^^ represents the Ricci scalar calculated using the three-dimensional metric hij and 
Eij is related to the extrinsic curvature of the constant time hypersurface and is given as 

1 



(16) 



E.. 



h,j - VjN, - ViNj 



(17) 



Varying the action (fT6il we get the constraint equation for and A^* as 



i?(3) -2V - N-^{EijE'^ -E'^)-N 



-2 i2 



0, 



(18) 

Ni + diifj where 



V, {Ei-SiE)]=0. 

Now we can decompose A^, into irrotational and incompressible parts as A^, 
(9jA^* = and expand A^, ■0 and A^* into powers of ( as 

N = 1 + ai + ^2 + 

= iv/'^ + iv/'^ + ... 

^ = + ^2 + (19) 

Using these expansions, the constraint equations f|T8|) can be solved order by order with 
metric f|T5|) and at first order one gets 



0, 



(20) 



Here d"^ = 5^Widj and the use of suitable choice of boundary conditions has been made to 

put A^^*^^'' = . As mentioned in 1^, to calculate the action up to n^^ order in we need 
to calculate A^ and Ni only up to the order-C"~^ and here the terms of order-^^ also drop 
out from the third order action, so equation fl20|) is sufficient to compute the action up to 
third order. So, after putting these solutions in ( fT6|l we get the action for second and third 
order in C as |lO, H Q 



^2 



dtd X 



dtd-^x 



(21) 



(22) 
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3.3. Two-point correlation function and power spectrum 

Now to calculate the two-point correlation function the quadratic part (!2Tll of the action 
is considered, which in conformal time {dr = ^) can be written as 

S^ = J drd^xah [C" - (5C)'] (23) 

Here / denotes derivative w.r.t conformal time r. The above action looks like an action of a 
massless scalar field in conformal spacetime and ( can be considered as the scalar field for 
quantization. ( can be written in terms of creation and annihilation operator as 

/d-^k ^ - /" d^k / ^ ^ X \ - 

^C(^, r)e^'-' = Jj^3 r)a^, + u\-k, T)al^) e*^^ (24) 

The equation of motion for C can be obtained by varying the action ( 123|) and is given by 



z 



C + 2-C - dX = 0. (25) 



z 



Here z"^ = 2a^e and we can define vj^ = z({k, r) and use equation (l24l) to get 



z" 



--]n = o- (26) 



The solution for the mode functions vj: can be obtained assuming Bunch Davies initial 
conditions and is given as 

= i fl - (27) 



/2A^ V 

Hence the basis function u{k, r) is 

uik, t) = ^ = ^= (1 + tkr) e-^'\ (28) 

The two point correlation function of the field Q in position space can be expressed as 

/d^kd^k' ^ ^ /- - N 

^^(0|C(fc,r)C(fc',r)|0)e<^'-+'='-^) (29) 

"^'^ u{k,T)\'e^^<^-^, (30) 



(2^) 



where we have used the relation {Q\C,{k,T)C{k' ,t)\Q) = {27v)'^6^ yk + k'j u{k,T)u*{—k' ,t). 
The power spectrum for C is defined by 

(0|C(^,r)C(fc"'',r)|0) = i2nfS' (k + k') P^{k) (31) 
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P^{k) = \u{k,T)\ (32) 
The another convention for the power spectrum, that is commonly used for data analysis, is 

A^ = |^KA;,r)| (33) 

In this case represents the variance of the classical fluctuations and the two-point corre- 
lation in position space becomes 

(C(f,r)C(f,r)) = j '^^yH^-y). (34) 

The power spectrum is calculated on super-horizon limit i.e. —kr « 1 in which = 
'^k Q-^kr gg|- ^Yie power spectrum as 

= (35) 

Now due to noncommutativity of spacetime the two point correlation function for field C 
gets modified [l2|. We will denote the field in noncommutative spacetime with a subscript 
6. Since here ( represents our quantum field, hence similar to (jl]) the twisted quantum field 
(0 can be expressed in terms of the untwisted field ( as 

(g(f,t) = ^(f,t)e5^MAP._ (3g) 

With the twisted quantum field one can compute the two-point correlation function in 
position space as 

{Ce{x,tKg{y,t')) = (C(f,t)e^^^.^^''C(y,t')e^^^''^'''') 

= (C(f,t)C(y,0)e-^^^-^^-, (37) 

where we have used the commutation relations between the field and the momentum operator 
[Pfj_, C] = —'id^(^. Now taking the Fourier transform on the right hand side we get 



(Pkd^k 
(27r)6 

'^-^^{Q\^{k^t)ak'MQ)e-^-^'^^^^^^^ 
|^(0|C(M)C(^,O|0)e(^^^^^'^^--^')e<^-^-^ 



(38) 
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Here & 
as 



. So the two-point correlation function in momentum space can be expressed 

(39) 

Now since in de Sitter space 



{0|C(fc,<)C(*',i')|0) = ei'-^'-iOIC I + ^1 C fk'.t' - ^ I |0> 



So in conformal time and in the hmit t' ^ t 

9^ -k' 



C\k,t' 



-k 



^C(k,re--^^ 



Hence the two-point function will be 



{ak,T)ak',T)) = (0|C(^fc,re-^— )C(fc',re^— )|0)e^ 

2 



kAk' 



u ( k, re^ 2 



{2TTf6^{k+k') 



Now we take the self-adjoint part of two-point correlation function defined as l3 



{0\Ce{k,r)Ce{k',r)\0)^ = - ( (0|G(A;, r)G(fc', r)|0) + (0|Ce(-A;, t)C.(-A;', r 



So the power spectrum can be obtained from ( l3Til as 

2 



P(eik) 



Since on super-horizon limit 



u ( fc, re^^ 2 



M — K, re 2 



1 / -i„-H 



V2k \ kr 



flO.fc 



6 2 SO 



P^^(A;) = P^(A;) cosh [He^-k 



(40) 

(41) 
(42) 



(43) 



(44) 



(45) 



(46) 



This power spectrum was derived in [12I, 1^| and was tested with the CMB observations in 
13| and it was found that the extra ^-dependent factor in the power spectrum can lead to 



direction dependent CMB anisotropy and future data at large / can tightly constraint the 
scale of noncommutativity. 
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3.4. Three-point function 



The primordial non-gaussianity in CMB arises due to the non-zero three-point and four- 
point correlation functions of curvature perturbations. These correlation functions were 
calculated for noncommutative spacetime in [l6| , where they have used SN formalism which 
ignores modifications to the correlation functions at Hubble crossing and also interaction 
between quantum fluctuations on sub-hubble scales with the super-hubble scale fluctuations 
at non-linear label. 

The third order action f l22|) obtained for C using ADM formalism is 



dt(rx 



(47) 



We put the value of i/j from Eq. fl2Ul) in this action, integrate by parts and use background 
Priedmann equations to get terms proportional to 



dtd-^x 



(48) 



where J-" = (r^C^ -|- terms with derivarives of () and ^ represents the terms propostional 
to the Gaussian action 5*2. We can again integrate by parts the above action to remove the 
terms involving dx and use the Gaussian field equation (l25l) to get 



5.3 



dtd X 



2 



(49) 



Now J-" = (?7 — e) + 2e9~^ (C^^C) and (9~^ is the inverse of and we have ignored the 
terms containing the derivatives of C in as they are negligible on super-horizon scales. 
One can get rid of the second term in the above action following field redefinition 



C ^ Cn + jC' 

After this field redefinition the three-point function becomes 



(50) 



(C(Xi)C(x2)C(a;3)) = {C,n{Xl)C,n{x2)C.n{xz)) + 



4 



{{C,n{xi)C,n{x2)){C,n{xi)Cn{xz)) + permutatious) 



+ 2^x^ {{C{xi)C{x2))dl^{Cn{xi)Cn{x3)) + peruiutations) . 

The first term in above expression represents the three point function while the last two 
terms represents the corrections to the three point function due to field redefinition. We 



(51) 
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will omit the subscript n in the following calculations. Now the interaction Hamiltonian to 
calculate the three-point function can be obtained from action ( 149|) i.e 



(52) 



As mentioned earlier, we use C, as the quantum field to compute the various correlation 
functions. Hence, to see the effects of noncommutative geometry on three point correlation 
function we replace the usual quantum field C, with the twisted field C,g both in the interaction 
Hamiltonian and in (|5T|) . Since the product of the twisted fields at the same spacetime point 
is given by the star-product [llj, the interaction Hamiltonian will be given as 



n{t') 



-2/ 



(53) 



Using the relation ^ between the twisted and untwisted quantum field and expression for 
star product ([5]), the interaction Hamiltonian becomes 



H(t') 



(54) 



The first term in the RHS of (ISTjl is computed using the in-in formalism |23| and is given by 



{Ce{xi)Ux2)U^,)) = -I / (it'(0|[G(xi)Ce(x2)Ce(x3),?/(t')]|0) 



to 



= -% \ dt'((0|O(a;i)O(a;2)O(x3)H(t')|0)-(0|H(t')O(a;i)O(x2)O(x3)|0)) 
Jto 

(55) 

Here the three point function is calculated at equal time i.e = ^2 = ^3 = ^- Initially we 
will write them differently for simplification but will put them equal before integration w.r.t 
t'. Let us now consider the first term of above equation with (|5^ and name it as (a). So 



a)=Ue' / dHa^E \ d'' x{^\Ce{xrXe{x2Xe{x^) C d~ \ 



62 



t'x 



|0) 



(56) 



Using the relation (1361) we can replace the twisted quantum fields in terms of the untwisted 
quantum fields and it gives 



(a) = Aie^ j dt'a^H j d3a;(0|C(a;i)C(2:2)C(a;3)e~^(^''^+^''^+^''^) 



X 62 



e2> 



t' ,X 



|0) 



(57) 



Aie^ j dt'a^H j rf3x(0|C(xi)C(2:2)C(a;3)e~^(^''^+^''^+^''^) 

i(^+£;+^)A^^g-i(^+^+£;)A^^g-f(^+£;+£;)A9^^-2^lQ^ 



X 6 2 



(58) 
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The above equation in Fourier space becomes 

1 fci(2vr)i8 



Here 



) = -Aie^ I dfa^'H [d^x [fl ^!^e<fci-^i+fe-x2+fc3-x3) 
X (0|c||^i,ti + 

X (o|c(^i,ti + 



i;i , ti H 







2 










h ^ ■ - 




■ h + -h 






2 










h ^ ■ ^4 - 







X (0|C A;i,ti + 



X C M^4,t - 



X c h,t' 







2 




















2 










+ ^■4 




■4 



X C yh, t' j |o) 

P = A 4 + A^2 A 4 + A^i A 4 + (a^i + k2 + (^4 + 4 + 4)) (60) 

Since we will express the three-point correlation function in momentum space, we can take 
the Fourier transform on both side of Eq. (1551) and take the limit ti = t2 = = t to get 

(a) = -z^*dt'((0|G(fci,t)O(^2,t)O(4,t)H(U')) 

= _4,e2y dt'a'H J d'x J ^-^^^{0\C{k,,t,)ah,t2K^^^ 
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where 



ti — t -\- 





2 














h + 






2 

^2 + ?- 




h + o^ 









2 




h + e^ -h + e^ 


■h 




2 






h + e^ -h + o^ 


■h 




2 






fci + ^ ■ 4 + ^ 


■h 



t2 = t + 
t3 = t + 
U = t' - 

h = t' - 

t, = t'- '"^^"'^ ^"^^^'^ '"^ (62) 

(63) 

A detailed calculation of this term is presented in Appendix and it is given as flA.lip 

(a) = t{2Ti)H^ (ki + ^2 + 4) n ^ K ^ ^''''"'■^ 

(64) 

Here K = ki + k2 + k^. Now similar calculations can be done for the second term in the 
three-point function f l55l) . Let us represent it as (6), 

[b) = I f dt'{0\n{t')Cg{x,)Ce{x2)Ce{x3m (65) 

J to 

The contribution due to this term in momentum space is given by flA.lGP 

, X ff^ 3 . £££^l%t£2±M l(fciAfc2+fc2Afc3+fclAfc3) 

(6) = e{27rfS' (k, + ^2 + 4) J] 0' k + ^^'^ 

So the contribution to the three-point function of C due to the first term of (15T!) in Fourier 
space is given as 



{Ce{h,t)Ce{k2,t)Ce{h,t)) = 2e{27T)'5' {k, + h + h) ^U^, 

i=l * 



COGh ^-^^-(^l+^^+fcs) ^ (fciAfc2+fc2Afe3+felAfc3) 

X — (^1^2 + perm.) 

(66) 

This concludes the calculations of the three-point function of the redefined field (n- Now 
to get the final three-point function of the field ( we need to consider the second and third 
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term of equation f l5T]) coming due to field redefinitions. The contribution to the three-point 
function due to first of these terms can be obtained using Wick's theorem and is given as 

(Ce(a;i)Ce(a;2)Ce(a;3) = ^^-^ ((Ce(a;i)Ce(a^2))(C0(a;i)Ce(2;3)) + perm.) (67) 

Now 



So 



//■ \r (rr W/r (t V (t W - I '^^^2^^fc3 1 f/^.ffca+fc3Liffc2+fc3V^l-»fc2-^2-»fc3-^3 

d^kid^k2d^k3 /-^ r \ H'^ 1 



(2-)7 (2^)9 [h + k2 + k,) 



16e2 A;|A;3 



^ ^He^-ki^-iki-xi~ik2-X2-ik3-X3 {^^) 

Here in the second step we have introduced a 6 function with integral over ki so that it 
matches with the results of the rest of the terms. So the contribution to the three-point 
function due to first field redefinition term in momentum space will be 

i=l * \ j / 

Now consider the second field redefinition term in Eq. ( ISTj) the contribution due to that is 
given as 

(C^(xi)Ce(x2)C<?(x3) = ^ {d-^' {{C9ix,)Ceix2))dl{Ceix,)Ce{x3))) + perm.) (71) 

Now 

dlACeixMx^)) = -j H^^^"^^-'^^^'^-^'^-'^^ (72) 

So 

{(:e[Xi}Q{x2}}d^^{(:e{xi}Q{xs}} - - J (27r)6 IQe^ k^k^ 

3 f d^kid^k2d^k3.3 . r . r \ 1 



^ ^H0^-ki ^-iki-xi~ik2-X2-ik3-X3 (73) 

So 

d^kid^k2d^k3 /-^ r\ 1 



5.-^ ((c.(xi)c.(x2))9,\(c.(xi)c.(x3))) = {2nr / fj; (^1 + ^2 + 4) 



^ ^He^^-ki^-iki-xi-ik2-X2-ik3-X3 {7"^) 
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So the contribution due to this term in Fourier space will be 



{Ce{h, t)Ce{h, t)Ce{h, t)) = ^i2nf5' (^k, + k, + k,^ J£^f[l^ [j^ k^k 



i=l 



(75) 



Now combining all the results from fl66|) . fl70|) . f l75|) for the various contributions to the 
three-point function of we get the final three-point function using Eq. flSTj) in momentum 
space as 

^4 3 



Where 



cosh ^-^^ + ^ ^«(feiAfc2+fc2Afc3 + felAfc3) 

^ = 4e ^ 



16e2 J- J. A;? 



(76) 



is: 



3^ 



(77) 



This is the main result of this paper. In the limit 6'^'^ — )■ the above expression becomes 
similar to expression for the three-point function in commutative spacetime (Eq. (4.5) and 
(4.6) of Maldacena lOj). Now due to translational invariance ki + k2 + = and on 
comparing our results with the commutative case [lO|, ll8| we see that the first term in (!77|1 
is modified due to a phase factor that depends on 6ij, while the second and the last terms 
are modified by exponential factors. The three-point correlation function here is complex so 
to see its observational effects we again take its self adjoint given as 16 



{Ce{ki,t)Ce{k2,t)Ce{ks,t)) 



M 



{Ceiki,t)Ceik2,t)C9{k3,t)) + {Cei-kut)Cei-k2,t)Ce{-k3,t)) 



+ 



{27tY5-' {ki + k2 + k3 
Tj — e 



3 I 



16e2 J- J- A:? 

i=l 



4ecos (^^^ 
K 



5Z 



J2 cosh (^h'^ ■ ^0 j + ^ {Y1 '^^^j (^^ ■ 



(7J 



4. IMPLICATIONS FOR OBSERVATIONS 



The non-gaussianity in CMB is described in terms of the angular three-point correlation 
functions in harmonic space called as "angular bispectrum", which is related to the three- 



dimensional bispectrum of the primordial curvature perturbations defined as |24j . l25l 



(C(A;i, t)C{k2, mk,, t)) = {2ny6' [k^ + k2 + k,]B^ {k,, k^, k,) 



(79) 
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We can generalize the above definition of bispectrum for the twisted quantum fields in 
noncommutative space time and it can be expressed using ( !78l) as 



3 I 



4ecos 



(fc)lAfc2 



K 



V 



J2 k- cosh (^h'? ■ 



+ 



hk] cosh [h^* ■ k 



Here the bispectrum also breaks the statistical isotropy. The anisotropic bispectrum also 
arises in the cases where the vector fields are also present during inflation jlsj. In [1^ 
the method to analyze these models in the light of new CMB data is derived. Current 
observational limits on non-gaussianity are given in terms of a non-linearity parameter /tvl 
that determines the amplitude and scale dependence of non-gaussianity. We define /at^ in a 
similar way as where it is assumed that the corrections to the standard power spectrum 
due to statistical anisotropy are very small. So 



5 



B 



NL 



6 P^ik^)P^,ik2) + P^ik2)Pciks) + Pdki)Pdh 
Using the power spectrum (1351) it becomes 

fclAfc2 



(81) 



NL 



6 kf 



cos 



4e- 



i<j \ i 



K 



+ ^ {^hk] cosh [h'? -k^ 



(82) 



k 



This kind of Jnl generally arises where the curvature perturbation is expressed as Q 

Cg + \Q cind Jnl peaks at the so called squeezed triangle limit defined as \ki\ = |^2| = 

and {k^l << k. So in a similar fashion the f^i for noncommutative case in the above limit 
is given as 



2e cos 



^ ^ - I (cosh [h^ . h) + cosh [h^ . k2 



(83) 



It is clear from the above expression that the amplitude of f^L is very small and of the 
order of slow-roll parameters for the case of small statistical anisotropy. But it has scale de- 
pendence and direction dependence that can help us to distinguish it from the commutative 
case. The current limits on the amplitude of Jnl are j23| Jni = 48 ±20 for squeezed triangle 
limit at 68% confidence level and the Planck is supposed to see the /nl ~ 5. The running 

of the non-gaussianity for ^ = in our case is similar to [16], and it is shown in 28. 29 
that we will be able to constraint n^G with a 1 — cr uncertainty of An^vG ~ 0.1 with ongoing 
Planck and future CMB and large scale structure observations. Hence even though the am- 
plitude of Jnl with the noncommutative geometry is of the order of slow-roll parameters, 
there is a possibility to see the effects of scale dependence of fm due to noncommutativity 
in future observations. 
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5. CONCLUSIONS 



Detection of primordial non-gaussianity in the CMB anisotropy and large scale structure 
is the main challenge of current and future observations and it can play an important role in 
discriminating various models of inflation. In this paper we have calculated the primordial 
non-gaussianity in single field inflation with spacetime noncommutativity. We have used 
Maldacena's approach [10] to compute the two-point and three-point correlation functions 
for the comoving curvature perturbation ( for the noncommutative case described by (l2| . 
Both the power spectrum and the bispectrum for this model are direction dependent and 
breaks the statistical isotropy due to the preferred direction of 6. This direction dependent 
power spectrum was analyzed by [l3| to put constraints on the scale of noncommutativity 
in the light of WMAP5 data and it was concluded that the WMAP5 data at high / is not 
sufficient to constraint the noncommutative scale 6. Since Planck will measure the CMB 
anisotropy at high-precession up to / > 2500, the author and collaborators plan to test the 
power spectrum fl?Il) after the release of Planck data. The statistical anisotropic bispectrum 
can be extracted from the three-point correlation function of CMB 2^ and for /jv^ ^ 30, 



future experiments could be sensitive to a ratio of the anisotropic to the isotropic amplitudes 
of the bispectrum up to 10%. The amplitude of the non-linearity parameter f^L for our case 
is very small for small statistical anisotropy but it has a high scale dependence. Ongoing 
Planck and future CMB and large scale structure observations would be able to measure 
the running of non-gaussianity up-to 1 — a uncertainty of Auj^g ~ 0.1 [29|, that can help us 
in distinguishing the effects due to noncommutativity. With the help of these observational 
signatures of noncommutative inflation we will, in future, be able to determine the scale of 
noncommut at i vity. 
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Appendix 

The integral appearing in the calculation of the first term in the three-point function (15 5 p 
in Fourier space can be read from Eq. f l61l) as 

(a) = -Ate^ J dt'a^'H j d'x J ^-^^L^{0\C{k,,t,)C{h^ 



X e 



where fjS are given by Eq. (1621) . Now we will calculate the six-point function entering in the 
above integrand separately and denote it as A. So 

A = (o|C(^i, ti)C(4, h)C{h. h)ah, U)C{h, hXih, t^M (A.2) 
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Now using Wick's theorem and leaving the disconnected diagrams we will get 6 terms in 
above expression. Let us consider one of them and denote it by Ai so 



A = (o| [e{h,ti),C{h,U)\ [e{k2,h),C{h,h)\ [e{h,t3),C{h,te)\ |o) (A.3) 

where the and denote the positive and negative frequency part of the quantum field 
C (see ([MD). Now since 

C+{h,ti),C{h,U)\ = {2n)''6^ (^h + h)u(ki,h^u' (A.4) 

we have 

Ai = {271^6^ + ^4) [h + h) [h + 4) u (^fci, ti) u (k2, ta) u (4, ^3) 

X ii* u) ii* (-h,t5) (-h, to) (A.5) 



Putting this back to the integral ( l59l) and denoting the contribution due to this term as (a)i 
and doing the delta integrals we get 

(a)i = -4ie^^ dt'a^H J d^x-^u (^i,t-^ u (^2,t-2^ u (^-^.t^ ii'' {ki,t^ ii^ 



-4:it\2'KY5'' [ki + k2 + h] / dt'a''H--^u[ki,ti]u[k2,t2]u[h,h 

Kg 



to 



-Pi 

3 2 ' 1 



(A.6) 



where Vi = k5=-k2 ke=-k3 ^^"^ calculated using these values of momenta. 

Here the limit of integration goes from to = —00 io t = 00. To solve the integral we will go 
to conformal time where the above limits correspond to r — > (— cxd,0). Now from Eq. (!28|) 
we have 



u{k, t) 



iH 



[1 + ikr) e 



-ikr 



(A.7) 



z y/Aek^ 

Since from Eq. f HOl) we know that the conformal time corresponding to tjS will be r x 
dependent term from Eq. (!62l) wc have terms like {t + 9 dependent term) for ti, t2, so for 
t — 7- 00 or r — i- conformal time corresponding to ti, t2, t^ will be zero. So in conformal time 
u (ki,ti] u ( ^2, ^2 ) u ( ^3, ^3 ) u \ki,0] u ( A;2, I u ( A^s, ) . Now we denote the conformal 



time corresponding to t' by r so ii* [kuU] ^ ''"*Vr and from (gO]) and we get 



T4 = T5 = tq = re" 



He'^ ■(ki+k2+k3) 



Now 



du* {ki,T4 



dr 



He^ ■(ki+k2 + k3) 



(A.8) 
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So the integral (\A.6\i becomes 



i=l * 



J— 00 



X 



1=1 * 

Here we have rotated the contour from (—00, 0) to i(oo, 0) and K = ki + k2 + k^. Now to 
calculate Vi let us recall ( BUj) 

V = kiAh + hAh + kiAh + kiA + 4 + 4) +4a (|a-4 + 4 + h^+hA + 4 + 4) 

Hence Vi = ki A k2 + k2 A k^ + ki A k^. Now rest of the terms in flA.2p can be found be 
different permutations of k^, k^, k^ and the phase factors will be same as Vi after imposing 
the different conditions due to delta function integrals. So from equation fl6Tl) we get the 
first term of the right hand side of equation fl55|) in Fourier space as 

(a) = e{2nf6^ [k^ + ^2 + 4) J] ^ k ^'''^^^ ^ ^^'''^'^ 

i=l * 

(A.ll) 

Now the second term in the three-point function (1551) is denoted as {h) and can be read from 
Eq. flXlal) as 



ft 

(b) = I I dt'{0\nit')Ce{x,)Ceix2)Ceixs)\0) 



to 

,2 



1 dt'a^H J 






j dt'a^H J 







t',x 



e 2 

t',x 



X C(a;i)C(a;2)C(2;3)|0)e-5(^^^+^^^+^^^) (A.12) 
Here we have used ( 136|) . Now in the Fourier space we get 



(b) = -Aie' [ dt'a'H [ d'x [ TT j!^ ^{kMk,ys^-itA{K+K+K) 

X (0|C(^4, t'Xih, t'nh, t')ah, tMh, t2)ah, h)\o) 

= Aie^ [ dt'a^'H [ d^x [ T\ ^!i_e<fc4+fc,+fe6)-V(^i-"^^+'^^-^'^+'^«-"^0 
J J -J i=i k&V''^) 

X (0|C(^4, t4)C(4, h)C,{h. tMh, hXih, t2)C(4, ta) |0)e5^ (A.13) 
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Here 



V = kiAk2 + k2Ak3 + kiAk3-(^ki + k2 + 4) A (k^ + h + 4) (A.14) 



0^ ■ {k2 + k^- k^- kr,- k( 



U = t + 



t2 = t + 



? ■ (^-ki + h-k4-h~h 
2 

^ ■ ^A;i — A;2 — ^4 — ^5 — fc( 



^-6 

ts = t + 



U = t' + 



h = t' + 



■{k^ + k2 + k 

2 

^ • (^1 + ^2 + k 



■{ki + k2 + h 



h = t' + ^ (A.15) 

Now all the calculations can be done for {h) as earlier and the final answer is 

(6) = e{2ixf5' (k, + ^2 + ^3) n 0- (^1^2 + perm.) 

i=l * 

(A.16) 
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